Monte Carlo simulations of magnetic ordering in the fee Kagome lattice 
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Monte Carlo simulation results are reported on magnetic ordering in ABC stacked Kagome layers 
with fee symmetry for both XY and Heisenberg models which include exchange interactions with the 
eight near-neighbors. Well known degeneracies of the 2D system persist in the 3D case and analysis 
of the numerical data provides strong evidence for a fluctuation-driven first-order transition to well- 
defined long-range order characterized as the layered q = (120-degree) spin structure. Effects 
of varying the inter-layer coupling are also examined. The results are relevant to understanding 
the role of frustration in IrMn3 alloys widely used by the magnetic storage industry as thin-films 
in the antiferromagnetic pinning layer in GMR and TMR spin valves. Despite the technological 
importance of this structure, it has not previously been noted that the magnetic Mn-ions of fee 
IrMn3 form Kagome layers. 

PACS numbers: 75.10.Hk, 75.30.kz, 75.50.Ee, 75.40. Cx 



I. INTRODUCTION 

The phenomenon of pinning the magnetization direc- 
tion of a thin-film ferromagnet due to coupling with an 
adjacent thin- film antiferromagent(AF) is key to spin- 
valve based magnetic transducer technology^ Under- 
standing exchange bias at the microscopic level has pro- 
gressed substantially in recent years and suggests that a 
mechanism to stabilize domains in the surface layers of 
the AF is essential. 3 Nearly degenerate energy states of 
the surface spin structure, such as found in geometrically 
frustrated AFs, can facilitate such domain formation. 4 
Although there are many features of the Ir-Mn com- 
pounds, such as high Neel temperatures (T/v ~ 900 K) 
which have contributed to their being one of the most 
commonly used materials for the pinning layer in mass 
produced spin valves, it has been argued that the spin 
frustration realized in the parent fee compound IrMn3 is 
relevant. 5 IrMn3 is in a class of magnetic compounds hav- 
ing the CuAu3 crystal structure which can be described 
as 2D triangular planes ABC stacked along (111) axes£ 
However, it has not been previously noted that the mag- 
netic Mn ions reside on sites in these planes that form the 
Kagome structure. The present work uses Monte Carlo 
(MC) simulations to explore the spin ordering and phase 
transitions with near-neighbor (NN) exchange interac- 
tions of such a 3D system, which we call the fee Kagome 
lattice. 

Heisenberg or XY spins on the 2D lattice formed from 
corner sharing triangles (see Fig. 1) exhibit a high de- 
gree of degeneracy in the case of near-neighbor exchange 
interactions where the only requirement is that the sum 
of spins on each triangle be zero (thus forming the 120°- 
spin structure) j2a£ The classical Heisenberg system shows 
entropy-driven planar spin nematic "order from disorder" 
at T = and has extensive entropy. The expected zero 



temperature limit of the specific heat, Cy = y| ks per 
spin, has been verified by Monte Carlo simulations^ 
Long range ground state Neel order has been demon- 
strated by adding further-neighbor exchange interactions 
in the Heisenberg case giving rise to so-called q = (see 
Fig. 1) or a/3 x a/3 spin structures^ The XY spin model 
on the 2D Kagome has received much less attention, 
likely due to the lack of relevant experimental systems. 
However, it has been argued that at T = the spin order 
should be that of the three-state Potts model, x i which is 
known to exhibit a weak first order transition at a non- 
zero temperature in 3D^ 




FIG. 1: (color online), (a) The fee Kagome lattice (solid cir- 
cles) formed by ABC stacking of triangular layers (solid and 
open cirlces) . In the case of IrMn3 , solid circles represent mag- 
netic Mn-ions and open circles denote non-magnetic Ir-ions. 
(b) Projection of the q — 0, 120° spin structure with ABC 
stacking along the (111) direction (out of the page) where 
colors distinguish the layers. 

Three-dimensional structures are formed from weakly 
coupled ABC stacked Kagome layers of magnetic ions 
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in the family of compounds with rhombohedral sym- 
metry known as the jarosites, AB3(S04)(OH)6, where 
a variety of experimental results suggest long range 
spin order of the q = 0, 120°-type, below tempera- 
tures in the range of 1 - 60 K j 13 i 14 Stacked Kagome 
layers have recently been investigated in an Fe-based 
metallo-organic compound which exhibits spin dynamics 
driven by frustration-induced domain walls. 4 Numerical 
simulations of the magnetic phase transitions in these 
systems have not been reported. Another 3D exam- 
ple is the hyperkagome spin lattice, derived from the 
pyrochlore structure by removing 1/4 of the magnetic 
sites, which has been shown to exhibit the degeneracy 
of the 2D Kagome lattice and MC simulations indicate 
a first order transition to octupolar ordering at very low 
temperature^ 

In addition to the very significant effort to understand 
the exchange bias mechanisms relevant to Ir-Mn thin 
films at the macroscopic level (see, e.g., Ref. [2] for a 
review), there have also been a number of works devoted 
to investigations of the antiferromagnetic spin order in 
bulk Ir-Mn (and sister) alloys at the microscopic level. 
Early neutron diffraction results on samples with vari- 
ous relative concentrations of Ir in the disordered form 
of I r Mni_ s alloys exhibit a simple two-sublattice mag- 
netic order which suggest a peak in Tn for values of 
x close to 0.25. 15 The magnetic structure of IrMn3 was 
found to possess the long-range 120°-type spin order be- 
low T/v — 960 K, referred to as the Tl state, which is 
equivalent to an ABC stacking of q = ordered Kagome 
planes. 16 This magnetic structure had previously been 
proposed for the ordered phases of RhMn3 and MnPt3 
alloys £ This conclusion is also consistent with first prin- 
ciples electronic structure calculations which, in addition, 
suggests that multiple-q spin stuctures are stabilized in 
the disordered alloys. 17 Most relevant to the present work 
involves a simulation of the spin structure of ordered 
IrMn3 using the stochastic Landau-Lifshitz- Gilbert equa- 
tions with a Heisenberg spin hamiltonian which includes 
a local (111) anisotropy term^ Such simulations are 
essentially equivalent to the MC approach. Model pa- 
rameters, which included a large number of exchange 
interactions up to 10 A between neighbors, were deter- 
mined from first-principles electronic structure calcula- 
tions. The q = (Tl) spin structure was verified and sim- 
ulations of the order parameter vs temperature yielded 
the estimate T/v — 1145 K. In the present work we 
examine the nature of the magnetic order, ground- state 
degeneracies and phase transitions in the fee Kagome 
Heisenberg and XY models with exchange interactions 
between the 8 NN shown in Fig. 2. This is achieved 
through extensive ground-state as well as Metropolis MC 
simulations. Results are presented for the temperature 
dependence of the internal energy, specific heat, order 
parameter and susceptibility all of which show clear in- 
dications of a phase transition to q = spin order for 
both models. Analysis of the energy histograms as well 
as finite-size scaling of the specific heat and Binder cu- 




FIG. 2: (color online). Illustration of the eight near- 
neighbors of the fee Kagome lattice projected onto the (111) 
plane (see Fig. 1). Four neighbors are in the (red) plane, and 
two each are in adjacent planes (black and blue). 

mulant indicate a strong first-order transition in the XY 
case but only weakly so for the Heisenberg model. We 
propose a model of the ground state degeneracy involv- 
ing planes of defects which appears to explain the mul- 
tiple values of the sub-lattice magnetization found from 
the simulations at low T. The model also explains the 
differences between MC simulation results performed as 
heating, cooling or independent temperature runs. Addi- 
tional MC simulation results are discussed for cases with 
weaker inter-layer exchange coupling which show the ex- 
pected decrease in T]y. The remainder of this paper is 
organized as follows. In Sec. II, the model hamiltonian 
is presented as well as details on how the simulations 
were performed. In Sec. Ill, the main results are shown 
for the temperature dependence of the various thermody- 
namic quantities. Finite-size scaling analysis results are 
also discussed. This is followed by a description of our 
model of degeneracies. Simulation results with weaker 
inter-layer coupling are presented in Sec. IV, followed by 
Sec. V where we discuss our results and directions for 
future simulations. 



II. MODEL AND SIMULATIONS 

Monte Carlo simulations were performed using the 
standard Metropolis algorithm onLxL Kagome planes, 
ABC stacked with L layers. Periodic boundary condi- 
tions were used on lattices with L=12, 18, 24, 30, 36, 
and 60. Between 10 5 and 10 7 MC steps (MCS) were 
used, with the initial ten percent discarded when cal- 
culating thermodynamic averages. Only NN exchange 
interactions were included in the simulations, as defined 
by the following Hamiltonian 

intra— plane inter— plane 

n = j Y. s, • s, • ./' ]T s, -s,, (i) 

i<j i<j 
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where the sum is over NN lattice sites, J = 1 represents 
AF coupling to the four in-plane NN's and J' > couples 
to the four NN sites in adjacent planes, as in Fig. 2. Most 
of the results below focus on the case of the fee lattice 
stucture where J' = J, giving the full eight NN interac- 
tions. Results are also presented for systems with weaker 
inter-plane interactions, J' < J. Simulations were per- 
formed assuming both XY and Heisenberg spin degrees 
of freedom. In the XY case, spins were assumed to lie in 
plane. 

Preliminary to performing simulations on the stacked 
Kagome lattice, extensive checks of the computer code 
were performed against published results for related sys- 
tems. The ground state spin configurations and energies 
for a 2D Kagome model which also included second and 
third neighbor exchange interactions reproduced results 
in Ref. [l0| for the boundaries between the q = and 
y/3 x y/3 phases. MC simulations on the 2D system for 
the specific heat at low temperature gave results consis- 
tent with those in Ref. 9 with NN exchange interations 
only. Additional simulations were made on the ordinary 
Heisenberg fee AF (ABC stacked triangular layers) where 
the temperature and first-order nature of the transition, 
based on the specific heat, were found to agree with pre- 
vious results. 19 



III. FCC KAGOME LATTICE 
A. Order of the transition 

A focus of these simulations was to establish the nature 
of the transitions to long range magnetic order in the XY 
and Hesienberg models on the fee Kagome lattice through 
evaluations of the energy and specific heat. Fig. 3 shows 
the energy vs temperature of these two models for the 
case of L = 24. Three types of simulation results are 
shown. In the case labelled cooling, the first simulation 
was run starting at a high temperature with a random 
initial spin configuration. Subsequent runs were made at 
decreasing temperatures where the initial spin configura- 
tion at each temperature was taken to be the final spin 
state of the previous run. A corresponding technique was 
used for heating runs, starting from the q = ground 
state. In these cases, thermal averaging was performed 
using 10 6 MCS on a single CPU. We also performed in- 
dependent temperature simulations where a new random 
initial state was assumed at each value of T. For this 
purpose, simulations at each temperature could each be 
performed on multiple, separate CPU and 10 7 MCS were 
used for averaging in these cases. The results show a 
clear indication of phase transition for both models and 
the discontinuity close to T/v = 0.760 for the XY model 
indicates it is first order for this two-component spin sys- 
tem. In the Heisenberg case, the order of the transition 
close to T/v = 0.476 is less clear and is further inves- 
tigated below. These results can be compared with the 
first order transition estimated to be at T/v = 0.446 in the 



ordinary fee Heisenberg AF which involves collinear long 
range magnetic order^ The fact that the three types of 
simulations result in identical energy curves will be rele- 
vant to the discussion below of degenerate spin states. 

Temperature 




FIG. 3: (color online). Energy of the fee Kagome (a) XY 
and (b) Heisenberg models for L = 24 with cooling, heating 
and independent temperature simulations. 

Energy correlations were used to calculate the specific 
heat for the two models, shown in Fig. 4, using the same 
simulation runs as for the energy. Well defined peaks oc- 
cur at temperatures corresponding to the features in the 
energy plots of Fig. 2 and again there is no difference in 
results obtained from the three types of simulations. The 
curves also appear to suggest a zero temperature limit for 
the specific heat of close to 1/2 k# per spin in the XY 
case and close to 1 for the Heisenberg model. The latter 
result is consistent with the value ^ ks per spin known 
to occur in the 2D Heisenberg system. 8,9 Such analysis 
would also suggest that the corresponding value for the 
XY model should be ^ per spin, again consistent 
with our results. Further indication that the degeneracies 
present in two dimensions persist for the 3D fee Kagome 
lattice are explored below. 

Additional simulations were performed to explore fi- 
nite size effects on the order of the phase transition for 
the Heisenberg case. An indication that this transition 
is weakly first order in the thermodynamic limit is sug- 
gested by the clearer discontinuity in the energy shown 
in Fig. 5 for the larger lattice L = 36 using 10 T MCS 
for averaging and independent temperature runs. Fur- 
ther evidence of the first order nature is also provided 
by the energy histograms shown in Fig. 6 which illus- 
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FIG. 4: (color online). Specific heat of the fee Kagome 
(a) XY and (b) Heisenberg models for L — 24 with cooling, 
heating and independent temperature simulations. 



trate a discontinuous jump in the energy minimum as 
the temperature varies only slightly around T/y for both 
the XY and Heisenberg models. 20 These histogram data 
also provide an accurate estimate of the respective transi- 
tion temperatures corresponding to the value of T which 
exhibits a double peak structure. Finally, the Binder en- 
ergy cumulant 20 was also calculated for the Heisenberg 
model using the larger lattice sizes with 10 T MCS for av- 
eraging, as shown in Fig. 7. The minimum exhibits clear 
finite-size volume scaling; however, its extrapolated value 
for L — > oo appears close to 2/3 expected of a continu- 
ous transition but our data are not accurate enough to 
make a conclusion on this point. It is noteworthy that 
for the frustrated stacked triangular AF, the extremely 
weak nature of its first order transition was confirmed 
by numerical simulations only after some twenty years of 
study by a large number of groups . 21 i 22 



B. Order parameter and susceptibility 

The q = ground state magnetic structure of the 2D 
Kagome lattice is defined by the spins on each triangle 
being 120° apart. This rule can also be satisfied in the 
ABC stacked fee structure with eight NN spins, as shown 
in Fig. 2. We have verified this to be the ground state 
for the 3D system and it is assumed here that this is the 
long range ordered state that occurs below the transition 
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FIG. 5: (color online). Energy of the fee Kagome Heisenberg 
model for L = 36 with independent temperature simulations. 



(a) 





T=0.755 

\ T=0.760 




\ T=0.765 























1.018 1.028 1.038 1.048 1.058 
- Energy 



(b) 



30000 
25000 
20000 





T=0.4760 

\ T=0.4763 




\ T=0.4766 / 































1.018 



1.038 
- Energy 



FIG. 6: (color online). Energy histograms for the fee Kagome 
(a) XY with L — 24 and (b) Heisenberg models with L=60. 



temperatures in the XY and Heisenberg cases considered 
in the previous section. The order parameter (OP) is 
defined through the three interpenetrating ferromagneti- 
cally aligned sublattices, (77 = 1,2,3) associated with 
the 120° spin structure 



M t = (3/7V){(M 2 + M 2 + Mi)/3}i 
where N is the number of sites and 



= (M*y + (My) 
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FIG. 7: (color online). (a) Binder energy cumulant of the 
fee Kagome Heisenberg model using independent tempera- 
ture simulations, (b) Scaling of the minimum value vs inverse 
volume. 



with 



Ml 



(4) 



where a = x,y,z. In the XY model, only x and y spin 
components are considered. 

Fig. 8 shows the temperature dependences of the cal- 
culated OPs for the XY and Heisenberg models with 
L = 24 for cooling, heating and independent temper- 
atures simulations using 10 6 MCS for averaging. The 
curves obtained by heating the system from its 3D q = 
ordered state show full saturation of the OPs for T — » 
and smooth, monotonically decreasing functions as T 
increases. The curves for the cooling process are also 
smooth but do not saturate at low T. The indepen- 
dent temperature runs are not monotonic and the OP 
value jumps between a number of values that mostly lie 
between the heating and cooling curves. These results, 
together with the fact there is no difference in energy 
between these three types of simulations shown in Fig. 3, 
indicate there are multiple degeneracies associated with 
the q = magnetic structure of the fee Kagome lattice. 
These points are explored further in the next section. 

Spin correlations were also studied through the OP 
susceptibility response function defined by 

X = ((Mi)-(\M t \f)/(k B T). (5) 

Fig. 9 shows cooling and heating runs only which are not 
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FIG. 8: (color online), q = order parameters for the (a) XY 
and (b) Heisenberg models with L = 24 from cooling, heating 
and independent temperature simulations. 



identical due to the degeneracies illustrated in Fig. 8. 
Independent temperature runs were also performed (not 
shown) and show large fluctuations due to the jumping 
between degenerate spin states (Fig. 8). The peaks in 
these curves occur near 0.76 and 0.47 for the XY and 
Heisenberg cases, respectively, are consistent with the 
transition temperature estimated based on the energy 
and specific heat anamolies, providing further evidence 
that the fee Kagome lattice indeed exhibits long range 
q = type spin order. 



C. Spin degeneracies 

In order to further investigate the impact of the q = 
spin configuration degeneracies associated with NN ex- 
change interactions on the stacked Kagome lattice, the 
three individual sub-lattice ferromagnetic OPs (Mi, M2, 
and M3) were calculated . Fig. 10 illustrates three exam- 
ple results of cooling runs performed on the XY model 
using three different random initial spin configurations. 
In each case, one of the OPs saturates at T = to the 
expected value of 1/3 and the other two OPs both tend 
toward the same smaller values, approximately given by 
0.22, 0.19, and 0.25. Additional simulations starting with 
different random initial spin configurations exhibit other 
saturation values of the OPs although none are less than 
about 0.16. Similar results were found for the Heisenberg 
model. 
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FIG. 9: (color online). Susceptibilities corresponding to the 
q — order parameter for the (a) XY and (b) Heisenberg 
models with L — 24 from cooling and heating simulations. 



These results can be understood by the observation 
that there is no change in energy with an interchange of 
the direction of two sub-lattice spins along a line in 2D or 
in a plane in 3D. This is equivalent to rotating the spins 
along a line (or plane) by 120°£^ The 2D case is illus- 
trated in Fig. 11. The possible values of the saturated 
sub-lattice magnetization due to this effect can be enu- 
merated by considering all possible defect lines/planes of 
sub-lattice switching, and is given by 



i^ 3 



(6) 



where n is the number of spins making the switch (e.g., 
between sub- lattices 1 and 2). In 2D the smallest devia- 
tion from full saturation of a sub-lattice occurs if there is 
only one row of switched spins so that n = L/2. In 3D, 
k planes of switching involves k (L/2) (L/2) spins up to a 
maximum where half of the population switches, n = |L 3 
(where the number of spins on each sub-lattice is \L 3 ). 
This means that the for the fee Kagome lattice the OP 
(at T = 0) of each sub-lattice will lie in the range [|, |] 
in the case of NN exchange interactions. The possible 
saturation values for the two switched sub-lattice OPs in 
the case L = 24 are shown in Fig. 12 and are consistent 
with the values found in the simulation results of Figs. 8 
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FIG. 10: (color online). The three ferromagnetic sub- lattice 
order paramaters associated with the q = spin state corre- 
sponding to three different cooling simulations, (a), (b), and 
(c), of the XY model with L = 24. 



and 10. 

Note that in each run, one of the sub-lattices under- 
goes no switching. It is possible that in different parts 
of the lattice different sub-lattices will take the role of 
the one that does not switch. Planes of switched spins 
formed in the cooling and heating processes usually are 
parallel; however, in some cases the planes may intersect 
each other. The spin configurations at these intersections 
impose an extra energy on the system. In a slow cooling 
process, these costly spin-configurations at the intersec- 
tions of the planes are eliminated by thermal relaxation, 
and finally all the switched planes are aligned parallel. It 
may be possible to trap these costly spin-configurations 
in a fast cooling process and values of sublattice magne- 
tization not predicted by Eq. (6) would occur. 



IV. INTER-LAYER COUPLING 

In addition to bulk materials which exhibit fee sym- 
metry, magnetic compounds composed of ABC stacked 
Kagome planes which are weakly exchange coupled are of 
interest (as noted in the Introduction). We examine here 
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FIG. 11: (color online). The three ferromagnetic sub- lattice 
spins shown as red, blue and black, (a) Perfect q = spin 
state, (b) Illustration of defects along the indicated horizontal 
line where two of the three ferromagnetic sub-lattice spins 
(black and red) are rotated by 120°. 



these curves are consistent with previous investigations 
of quasi-2D systems j 23 i 24 
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FIG. 12: (color online). Possible values of the sub-lattice 
order parameter due to switched spins in the case of L = 24. 



the dependence of the transition temperature of both 
XY and Heisenberg models on the inter-layer exchange 
coupling J' through simulations of the specific heat on 
L = 24 systems with periodic boundary conditions and 
using 10 7 MCS for averaging. Values 0.05 < J' < 0.95 
in steps of 0.05, with intra-plane exchange J = 1, were 
considered in independent temperature simulations dis- 
tributed over 400 CPUs. 

Example results for a few values of J' are shown in 
Fig. 13 and a summary of the dependence of T/y on J' 
is given in Fig. 14. The absence of a detectable peak 
in the specific heat at the smallest J' over this range of 
T > 0.05 is consistent with the predicted lack of long 
range spin order in 2D. We find that the specific heat 
tends to the T — > values observed for the fee case (Fig. 
4), independent of J' . This result again suggests that NN 
inter-layer coupling does not alter the nature of the spin 
degeneracies of the 2D system. The general shapes of 
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FIG. 13: (color online). Example specific heat simulation 
results for interlayer coupling J' <1 as indicated for the (a) 
XY and (b) Heisenberg models from independent temperature 
simulations with L = 24. 
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V. SUMMARY AND CONCLUSIONS 

A signifcant conclusion from the results presented in 
this work is that the fee Kagome spin lattice with NN 
exchange interactions (J) exhibits long range magnetic 
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order of the q = (120°) type through first order phase 
transitions at temperatures 0.760 J and 0.476 J for XY 
and Heisenberg models, respectively. In the Heisenberg 
case, finite-size scaling of the energy cumulant, along 
with other results, suggest that the first-order nature 
of the transition is very weak. For both of these mod- 
els, mean field theory would predict a continuous phase 
transition. Spin configuration degeneracies well known 
in the 2D system as lines of defects appear to persist 
in the 3D case as planes of defects and lead to multi- 
ple (but well defined at T = 0) possible values of the 
sub-lattice magnetization order parameters. These ob- 
servations suggest that the transitions are driven by the 
order-by-disorder phenomenon. Additional simulations 
of ABC stacked Kagome planes with weaker inter-layer 
coupling, J', exhibit similar types of transitions and de- 
generacies, with T/v decreasing monotomically to zero as 
J' 0. 

These results are relevant to a number of experimen- 
tal systems. Fee IrMn3 is one of the Ir-Mn compounds 
well known in the magnetic recording industry as useful 
for the pinning AF layer in spin valve devices. Previous 
neutron scattering experiments have confirmed that this 
and sister compounds show long-range q = spin or- 
der (known in the industry as Tl magnetic order). The 
pinning of the magnetization in a ferromagnetic layer 
through exchange coupling with an adjacent AF is be- 
lieved to be due to the presence of defects and domains. 
The results of the present work emphasize that geomet- 
rical frustration is useful in generating these spin defects 
in compounds with the fee Kagome spin structure, like 
IrMn 3 . 



This work represents a prelimnary investigation of a 
number of problems related to both fundamental aspects 
of geometrical frustation as well as exchange pinning. 
We have already performed initial simulations on the be- 
haviour of the fee Kagome system in an applied mag- 
netic field. Such studies previously revealed a tricritcal 
point in the case of the stacked triangular lattice and 
further illuminated the nature of the first order transi- 
tion for that lattice structure^ 2 Another extension of the 
present work is to examine the impact of the single-ion 
anisotropy noted in Ref. [18], as well as Dzyaloshinskii- 
Moriya interactions, 25 on the nature of the phase tran- 
sitions in 3D. MC simulations of thin film ABC stacked 
Kagome systems which include anisotropy as well as sur- 
face effects (on both exchange and anisotropy) would be 
useful to study surface spin-states and their possible im- 
pact on defects and domain formation. Extensions to 
include also a ferromagnetic layer with dipole interac- 
tions would reveal more about the microscopic mecha- 
nisms important for exchange pinning for In- 
vestigations of spin waves in stacked Kagome systems^ 
would compliment previous studies of the layered trian- 
gular AF 27 and would provide insight into the mechanism 
relevant to high frequency spin- valve sensor response. 
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